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Two components of the error of a troposphere calibration measurement have been 
quantified by theoretical csdculations. The first component is a beam mismatch er- 
ror, which occurs when the calibration instrument senses a conical volume different 
from the cylindrical volume sampled by a DSN antenna. The second component is a 
beam offset error , which occurs if the calibration instrument is not mounted on the 
axis of the DSN antenna. These two error sources were calculated for both delay 
(e.g., VLBI) and delay rate (e.g., Doppler) measurements. The beam mismatch 
error for both delay and delay rate drops rapidly as the beamwidth of the tropo- 
sphere calibration instrument (e.g., a water vapor radiometer or an infrared Fourier 
transform spectrometer) is reduced. At a 10-deg elevation angle , the instantaneous 
beam mismatch error is 1.0 mm for a 6-deg beamwidth and 0.09 mm for a 0.5-deg 
beam (these are the full angular widths of a circular beam with uniform gain out to 
a sharp cutoff). Time averaging for 60-100 sec will reduce these errors by factors 
of 1. 2-2.2. At a 20-deg elevation angle, the lower limit for current Doppler obser- 
vations, the beam-mismatch delay rate error is an Allan standard deviation over 
100 sec of 1.1 x 10 ~ 14 with a 4-deg beam and 1.3 x 10 ~ 15 for a 0.5-deg beam. 

A 50-m beam offset would result in a fairly modest (compared to other expected 
error sources) delay error (<0.3 mm for 60-sec integrations at any elevation an- 
gle >6 deg). However, the same offset would cause a large error in delay rate 
measurements (e.g., an Allan standard deviation of 1.2xl0~ 14 over 100 sec at a 
20-deg eievation angle), which would dominate over other known error sources if 
the beamwidth is 2 deg or smaller. An on-axis location is essentia/ for accurate 
troposphere calibration of delay rate measurements. A half-power beamwidth (for 
a beam with a tapered gain profile) of 1.2 deg or smaller is desired for calibration of 
all types of radio metrics. A water-vapor radiometer calibration beam of this size 
with very low sidelobes would require a clear aperture antenna with a diameter of 
at least 1.5 m if the primary water vapor sensing channel were in the 20-22 GHz 
range. 
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I. Introduction 

A. Background 

Spacecraft tracking with the DSN involves radio metric 
(delay or delay rate) measurements over the spacecraft- 
ground path, sometimes with similar measurements of nat- 
ural radio sources for purposes of calibration. Any non- 
geometric delay or delay rate will (if not perfectly cali- 
brated) degrade the accuracy of the tracking measurement. 
With the increase in the radio frequency used for space- 
craft tracking from 2.3 GHz (past missions) to 8.4 GHz 
(current missions) to 32 GHz (future missions), the effect 
of charged particles along the spacecraft-ground path is 
becoming much less significant (Ti on oc j/“ 2 , where r lon is 
the delay due to charged particles and v is the radio fre- 
quency of the link; the same proportionality holds for delay 
rates). The effect of the neutral (and nondispersive) tropo- 
sphere is consequently becoming more dominant. For 8.4- 
GHz single- frequency observations, the contributions from 
charged particles and from the troposphere are compara- 
ble. For 32-GHz single-frequency tracking (future space- 
craft missions) or for dual-frequency measurements (cur- 
rent very long baseline interferometry [VLBI] tracking), 
the troposphere will dominate over other media effects. 
The radio metric data themselves can be used to solve 
for the tropospheric delay. However, spatial and temporal 
irregularities in the troposphere limit the accuracy with 
which this process can be done. 

A VLBI tracking demonstration on a natural source [1] 
showed that the tropospheric delay in the direction toward 
a spacecraft could be estimated to an accuracy of approx- 
imately 1 cm, giving an angular measurement accuracy of 
1 nanoradian (nrad). Troposphere calibration is essential 
in order to improve the accuracy of radio metric tracking 
beyond this capability. For a goal of an order-of-magnitude 
improvement in VLBI tracking accuracy to 100 picoradi- 
ans (prad), the contribution of the troposphere must be 
reduced to 1 mm, through a combination of calibration 
and parameter estimation. 

Doppler (delay rate) spacecraft tracking measurements 
are used both for navigation and for gravitational wave 
searches. The troposphere calibration goals for future 
Doppler navigation using Ka-band (31-34 GHz) have not 
yet been completely specified. However, they are approxi- 
mately 1 x 10“ 15 in an Allan standard dev iation over 3000- 
10,000 sec. 1 For the Cassini Gravitational Wave Search 
experiment (the first one to use a Ka-band link), the goals 


1 S. Thurman, personal communication, Navigation Systems Sec- 
tion, Jet Propulsion Laboratory, Pasadena, California, February 3, 
1993. 


for troposphere calibration are Allan standard deviations 
of 4 x 10“ 16 at 1000 and 10,000 sec and 1.3 x 10 -15 at 
100 sec. The experiment requirements are a factor of 4 
less stringent. 2 The goals for this experiment were chosen 
so that the contribution to the delay rate from the uncali- 
brated troposphere would not dominate the total error at 
any time scale. 

In this article, the terms “delay” and “delay rate” refer 
to characteristics of the troposphere above a single DSN 
antenna. This is different than one common use of these 
terms in DSN spacecraft navigation, where “delay” and 
“delay rate” refer to the difference (or its time derivative) 
in arrival times at two DSN antennas (which are often 
widely spaced). For such station-differenced radio metric 
data types, the effects calculated here apply to both sta- 
tions. For widely spaced stations (e.g., at different DSN 
sites), the delay and delay rate effects in this article should 
be added in quadrature. 

B. Beam Mismatch and Offset 

Several types of remote-sensing troposphere calibration 
systems have been constructed or proposed. These include 
water vapor radiometers (WVR’s) [2], infrared Fourier 
transform spectrometers (FTS’s) [3], and tracking of global 
positioning satellites (GPS’s) [4]. For all these systems, 
the delay is measured (or inferred from thermal emission) 
over some volume of troposphere. This volume is ap- 
proximately conical for techniques that measure emission 
(WVH, FTS) and cylindrical for techniques that measure 
delay or absorption (GPS). In both cases, the measured 
volume of troposphere will not match the volume sampled 
by a DSN antenna, which is cylindrical with a diameter 
equal to the DSN antenna diameter, axis centered on the 
antenna, and extending in the direction of the spacecraft 
or natural source that is being observed. This cylindrical 
shape results from the circular shape of the DSN aperture 
and the fact that the observed radio signal arises from a 
source on the sky that has a very small (<Cl-arcsec) angu- 
lar extent. The volume mismatch can be broken down into 
two components: a mismatch in the beam shape and/or 
size (e.g,, conical versus cylindrical) and a horizontal off- 
set in the axes of the two instruments (resulting from the 
calibration instrument being mounted on the edge of the 
DSN antenna or off the antenna entirely). There are in- 
homogeneities in tropospheric refractivity, predominately 
due to variable water vapor density. As a result, both a 
beam mismatch and a beam offset would be expected to 


2 S. Dolinsky, “Cassini Ka-Band Radio Science Study,” (interned 
document), Jet Propulsion Laboratory, Pasadena, California, Au- 
gust 1992. 
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cause significant errors. The purpose of this article is to 
quantify these errors. 

Other error sources in the calibration process (e.g., 
the brightness temperature-to-path delay conversion for 
water-vapor radiometer measurements) will be important 
and are likely to be larger than the beam-mismatch and 
beam-offset errors in some cases. However, these error 
sources have been neglected in this article, so that the ef- 
fect of the beam- mismatch and beam-offset errors can be 
more fully studied. In any case, these other (e.g., non- 
beam-related) error sources could potentially be reduced 
to very small levels by advanced instruments and data- 
processing techniques. 

II. Definitions and Calculation Methods 

A. Basic Definitions 

The wet tropospheric delay tdsn(xdsn jO sampled at 
time t by a DSN antenna of diameter dosN whose vertex 
is at position Xdsn and which is pointed at elevation angle 
9 dsn and azimuth <j>DSN is 

tdsn(^dsn^) = — 72 / RdR / duj 

^^DSN Sin “OSN JO Jo 

X f dz X(XDSN 
J — Jtc OBWCOS &DSN 

+ r(0DSNi<l>DSN>R>ViZ)>t) ( 1 ) 

where is th e refractivity (x s n — 1, where n is 

the index of refraction) associated with water vapor at 
position x and time t ; v(0dsn, <f>DSN, #,u;, z) is the posi- 
tion, relative to the antenna vertex, of a vector in direc- 
tion (Odsn, </>dsn) from polar coordinates ( R,u > ) on the 
antenna to height z (the antenna vertex is defined as ly- 
ing at z = 0); and h is the height of the wet troposphere 
slab, above which x is assumed to be negligible. The DSN 
aperture has been approximated by a flat disk. In Eq. 
(1), the delay has been integrated across the DSN aper- 
ture. The actual delay measured by a DSN antenna is 
derived from the phase of a spatially integrated complex 
signal. However, because the signal is coherent across the 
DSN aperture (i.e., the phase variations are <<1 rad), the 
phase of the integrated signal is equal to the integral of the 
phase across the aperture. 


The wet troposphere delay r ca f(i? ca /, t) measured at 
time t by a calibration instrument with a conical beam 
at location x ca i is 

. i . , r du 

r C ai(x cai ,t) = — sin udu J r 

n Jo Jo Sin 0(u,w) 

X / dzx(xcai + r{e,<j>,z),t) (2) 

Jo 

The calibration-beam cross section is assumed to be coni- 
cal, with uniform gain out to a full angular width 0. The 
correspondence between the full width of such a beam and 
the half-power beamwidth (HPBW) of a Gaussian beam 
profile is discussed in Section III. A. The beam solid angle 
is Q = 27t[1 — cos (0/2)] and r{9 , 0, z) is the position at 
height z } relative to the calibration instrument, of a vector 
along elevation angle 9 and azimuth 0. 

In order for the difference between tdsn(zdsn ,t) and 
T C al(&cai,t) to have zero mean, the axis of the calibration 
beam must be pointed at an elevation angle 9 ca \ that is 
higher than Odsn • There is a nonlinear gradient of air 
mass A with elevation angle 9 across the calibration beam, 
so that the mean air mass corresponds to an elevation 
angle lower than 9 ca i. [The air mass A{6 ) is the tropo- 
spheric path length at elevation angle 0, normalized to 
the zenith path length, for a homogeneous troposphere. 
For the plane-parallel, zero-bending approximation used 
in these calculations, A(9) — l/sin#]. This effect has 
been previously noted [5]. For the circular calibration- 
beam cross section used here, 9 ca \ was calculated numer- 
ically. For the cases in this article, the largest value of 
9 ca i - ®dsn was 0.23 deg (for Odsn — 10 deg and 0 = 
6 deg). For a more typical case, the difference in elevation 
angles was much smaller (e.g., for 6 dsn = 20 deg, 0 — 
1 deg, 9 c<x j - 6 dsn — 0.003 deg). 

The absolute value of the average instantaneous differ- 
ence Ar <njt between the tropospheric delay sampled by 
the DSN antenna and that measured by the calibration 
instrument at time t is 


Ar,* n4 i(t) = yJ([TDSN(XDSN, t) - T ca t(xDSN +/?>0] 2 } 

(3) 

where p is the vector separation in the horizontal plane 
between the axis of the DSN antenna and the calibration 
instrument. When delays are averaged over a time At (as 
in a VLBI measurement), the delay difference A r avg is 
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AT av g(t 0 , A^) 


D x (x i,xi - x 2 ) = 


C]^_ 

1 + (R/L, at ) 213 


( 6 ) 


2 

I (4) 


where t 0 is the epoch at which the integration begins. For 
all cases in this article where Eqs. (3) and (4) were evalu- 
ated, 0 ca j — Odsn was chosen so as to give a zero mean for 
the delay difference. Therefore, Eqs. (3) and (4) yielded 
the standard deviation of the delay difference. 

The Allan standard deviation <r y (Af) [6] of the delay 
difference tpsn(%dsn >0 ~ T ea i(xDSN + P > 0 over a ^ me 
interval At is 


1 

At 




[At 

/ T D SN{XDSN,to + 0 

Jo 


dt 


r 

- / T ca i(x DS N + + 0 

Jo 


dt 


O-y(At) 



+ 2A*) - 2Ar, •„,«(< + At) + A W (t)] 2 ) 
2(A<) 2 


(5) 


0 - y (A*) will be used as a measure of the delay rate error 
due to beam mismatch and offset. 

B. Method of Calculation 

In order to evaluate Eqs. (3)-(5) for given values of 
dpsNy Odsn > <i>DSN> and the expressions were all ex- 
panded to integrals of x» using Eqs. (1) and (2), Prod- 
ucts of integrals were converted to multidimensional inte- 
grals, and the order of integration and expectation value 
was reversed. This resulted in multidimensional inte- 
grals of (x 0 ?i)x(^ 2 )}> where and x 2 are two differ- 
ent vector locations above the DSN antenna and/or the 
calibration instrument. These integrals were then con- 
verted to integrals of a structure function D x (x i, £i - z 2 ). 
(d x (x 1 , X\ - x 2 ) = (\x(x l) - x(*a)] S )) • This procedure 
has been described elsewhere [7,5]. 


where R = ||xi - x 2 \\ and L aa t is a saturation length 
needed to make D x converge at large distances. For the 
calculations reported here, the values of the parameters 
adopted were C x — 1.1 x 10“ 7 m” 1 ^ 3 , L $a t = 3000 km, 
v w — 8 m/sec, and the slab height h = 2 km. These val- 
ues give the best match to observed 200-sec delay rates 
and observed daily zenith delay scatter at DSN sites [1]. 
Connected-element radio interferometric observations over 
a 20 km-baseline (DSS 13-DSS 15) at Goldstone provide 
one test of the range over which C x can vary. 3 A series of 
10 such observations over a period of 2-3 years was ana- 
lyzed using the model of [7] to fit the observed scatter in 
residual delays as a function of time separation, angular 
separation, and elevation angle. One parameter adjust- 
ment of C x resulted in a \ 2 agreement of approximately 1 
per degree of freedom in all cases, with the best-fit value 
of C x always within a factor of 2 of the value used in this 
article. 

The multidimensional integrals of D x were evaluated 
numerically, using Simpson’s Rule. The integration step 
sizes were reduced until convergence was obtained. 

C. Validity of Tropospheric Model 

Evidence of the validity of this tropospheric model 
(Kolmogorov turbulence, frozen flow) comes from several 
arguments. The model prediction for the dependence of 
delay rate (at a 200-sec time scale) upon elevation an- 
gle agrees with the observed shape from 3 years of DSN 
VLBI data [7]. The phase structure function from non- 
DSN VLBI observations [8] shows the same power law 
dependence (with a curvature at integration times near 
200 sec) as that predicted by the model. The delay struc- 
ture function as measured by a solar hygrometer at JPL 4 
shows approximately the same power law dependence at 
short and long time scales as for the model. In addition 
to the connected-element radio interferometric data (on 
natural sources) cited above, additional observations of 
the Ulysses spacecraft with the same interferometer have 
been performed. 5 When C x is adjusted so as to fit point- 
point delay scatter of adjacent time points, the estimated 
difference in the zenith delays between the two antennas 


The model of [7] was used to evaluate these integrals. 
This model assumes that time variations in x are due t° 
spatial variations which move past the observer at a con- 
stant vector wind velocity v w . It assumes that the wet tro- 
posphere has constant average properties up to a height h, 
with zero refractivity above. D x has the analytical form 


3 C. Edwards, personal communication, Tracking Systems and Ap- 
plications Section, Jet Propulsion Laboratory, Pasadena, Califor- 
nia, May 28, 1993. 

4 R. Treuhaft, personal communication, Tracking Systems and Ap- 
plications Section, Jet Propulsion Laboratory, Pasadena, Califor- 
nia, March 1993. 

5 C. Edwards, op. cit. 



(DSS 13 and DSS 15) was reproduced within a x 2 of 1 per 
degree of freedom. 

The tropospheric slab approximation (constant struc- 
ture function from the surface up to a cutoff altitude) was 
tested numerically. As described in Section III. A, A r,* n , t 
was calculated for two {Odsn ><j>DSN , i>) cases with a more 
realistic structure function, one in which the refractivity 
variations decrease exponentially with height. In these two 
cases, ATi ns t differed from the slab model value by 3 and 
11 percents. 


III. Results: Delay Errors 

A. Beam Mismatch Effects 

This subsection presents the errors which occur when 
the calibration instrument is mounted on the axis of the 
DSN antenna. The errors are due to the mismatch be- 
tween the conical beam of the calibration instrument and 
the cylindrical beam (a good approximation for the tro- 
pospheric volume sampled when observing a point source) 
of the DSN antenna. For the case where the conical cali- 
bration beam, projected to a typical atmospheric height, 
subtends a much larger diameter than the DSN beam, this 
error is sometimes referred to as beam smearing. The in- 
stantaneous and time-averaged errors are given by Eqs. (3) 
and (4), respectively, with p = 0 in both cases. 

1. Zero Integration Time (Instantaneous Beam 
Mismatch Error). The simplest approximation used was 
a pencil DSN beam (i.e., a cylinder of diameter zero), a 
conical calibration beam, and zero integration time. Re- 
sults are shown in Fig. 1. Note that 9 cal - 9 DS n > 0 
in order to remove the systematic offset between r ca j and 

T DSN- 


Seven of the values shown in Fig. 1 were recalculated, 
using a 34-m-diam DSN antenna beam. The physical mis- 
match between such a beam and a conical calibration beam 
was smaller than for a pencil DSN beam. When the max- 
imum calibration-beam diameter {\ph/ sin Oosn) was not 
substantially larger than the DSN beam diameter ( h is the 
wet troposphere slab height), there was a significant reduc- 
tion in the delay calibration error At. (For these cases, 
the DSN cylindrical beam is a significantly better match 
to the conical calibration beam than is a pencil beam.) 
However, these were the cases where the delay error was 
small (<0.1 mm) with a pencil beam calculation. When a 
pencil beam calculation yielded a delay error >0.1 mm, the 
delay error with a cylindrical DSN beam was only slightly 
smaller than with a pencil DSN beam. 


For instantaneous (Af = 0) results, no transport of 
refractivity irregularities occurs. The results represent 
purely spatial irregularities and do not depend on the wind 
velocity or direction. 

2. Nonzero Integration Time. For four cases, the 
beam mismatch error was calculated for integration times 
of 10 and 100 sec, with a 34-m-diam DSN beam. For 
nonzero integration times, some of the refractivity irreg- 
ularities move through both beams, reducing the beam 
mismatch error. A wind velocity of 8 m/sec was assumed, 
with an azimuth 45 deg different from the DSN pointing 
direction (the dependence of the time-averaged error upon 
wind direction was not tested, due to limited computer 
time). The results are shown in Fig. 2. Delay integration 
for 100 sec reduces the beam mismatch error by factors 
of 1.2-2. 2. The smaller factors correspond to cases where 
the instantaneous beam mismatch error was largest. For 
these cases, the DSN beam cross-sectional area at most 
heights is much smaller than that of the calibration beam. 
Therefore, most of the refractivity irregularities in the cali- 
bration beam which are transported by the wind are never 
seen by the DSN beam, and time averaging has little effect. 

For an overall line-of-sight calibration accuracy goal of 
1 mm, it is desirable that any one error component con- 
tribute no more than 20 percent to a 1-mm 2 total variance. 
T his y ields a maximum allowed beam-mismatch error of 
\/0.20 mm = 0.45 mm. For a 100-sec integration time, a 
1-deg full-width calibration beam is acceptable at all ele- 
vation angles, and a 2-deg full-width beam is acceptable 
for elevation angles larger than approximately 7-8 deg. 

3. Test of Wet Troposphere Slab Approxima- 
tion. As a test of the validity of the wet troposphere slab 
model, the instantaneous delay error with a pencil DSN 
beam was recalculated using a modified structure function 
for two sets of (OdsNj^P) parameters. This structure func- 
tion was adopted from previous work 6 and had the form 


D x (x u x x - x 2 ) = 


C 2 ft 2 / 3 

1 + (R/L sat )V* 


e — (*i +*a)/ *o 


( 7 ) 


where ft = |[£i — ^2 1 1 f ^1 and z 2 are the heights of the two 
points Xi and x 2i L sat is a saturation scale length, and z 0 
is a wet troposphere scale height. This structure function 
accounts for the falloff of water vapor density with height. 


6 J. J. Wiseman, “Altitude Dependence and Saturation Modifica- 
tions to the Treuhaft-Lanyi Model of Tropospheric Delay Fluctua- 
tions,” JPL Interoffice Memorandum 335.3-87-203 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, December 
21, 1987. 
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Values of C x = 1.34 x 10“ 7 m- 1/3 , L sat = 3000 km, and 
zq = 2 km were used. 7 The upper height of integration 
was 10 km. The results (for zero integration time and a 
pencil, or zero-diameter, DSN beam) were 0.15 mm for 
OdsN = 10 deg and ij? = 1 deg, and 0.066 mm for 6 dsn 
— 30 deg and ifi = 2 deg. These results differed from the 
corresponding values in a slab model by 11 and 3 percents, 
respectively. This level of agreement suggests that the use 
of a wet troposphere slab approximation was a reasonable 
one for the purpose of calculating beam mismatch errors. 

4. Application to Beams With Tapered Gain 
Profiles. The beam mismatch results presented in Figs. 1 
and 2 are for a uniform calibration beam gain out to a 
sharp cutoff at rp/2. Actual beam profiles will be tapered, 
with a Gaussian profile often used for analytical calcula- 
tions. A Gaussian profile was not used for the calculations 
in this article because it would have required an excessive 
amount of computer time. 

An approximate correspondence between Gaussian 
HPBW and the full width ^ofa uniform gain beam pro- 
file used in these calculations can be made. In Fig. 1, 
the dependence of the beam mismatch error Ar* n ,* upon 
beamwidth is Ar injt oc xp exp , where exp «1 at low 
elevation angles {9 dsn = 6 or 10 deg) where the error 
is of most concern (i.e., largest). For a circular beam 
with a uniform gain out to a sharp cutoff, the mean 
value £ uniform for the angular separation from the axis 
is ^uniform = $ f 3. For a circular beam with a Gaus- 
sian profile, the mean value ^Gau ssian far this separa- 
tion is ^Gaussian = HPBWv/i/4 M 0.532 HPBW. 
Assume that Ary na * oc ip exp =>• Ar,- n5 * oc £ exp . For 
A Ti n8 t(ttniform) = A T inst (Gaussian), {uniform = 

^Gaussian , which yields HPBW = 0.63^. Therefore, beam 
mismatch errors in Figs. 1 and 2 will apply approximately 
to HPBW = 0.63^>. A roughly similar relation will hold 
for beam-mismatch delay rates (Tables 1-2, discussed in 
Section IV. B). For both delay and delay rate, the corre- 
spondence factor will differ very little from 0.63 at higher 
elevation angles (e.g., the factor is 0.64 if A Ti nst oc t/> 0 8 , 
as observed at 6 DS n = 00 deg for delay errors). 

5. Tests of Numerical Computations. The ac- 
curacy of the numerical integrations used to generate the 
beam mismatch results was tested in two ways. The first 
test used a quadratic structure function. 

D X (R) = C\R> (8) 


7 Ibid. 


Using this structure function, the integrals could be eval- 
uated analytically in many cases 8 (e.g., with pencil DSN 
or calibration beams, cylindrical DSN beams, or fan cal- 
ibration beams: a nonzero angular extent in only one di- 
mension). The analytical and numerical beam mismatch 
results for the quadratic structure function agreed to seven 
significant figures. 

The second test consisted of independent numerical in- 
tegrations, using the structure function of Eq. (6). These 
calculations used a code written for the formalism in [5], 
corrected for a numerical error [9]. A diverging beam with 
a square cross section and uniform gain was used in these 
calculations. For a square beam, £ uniform(square beam) = 

0.383^ = lAbi U niform(circular beamy Given the as- 
sumed relationship between Ar,* n ,< and £ discussed above, 
one would expect A Ti nst (square beam) « 1.15Ar* n ,t 
(circular beam) for the same value of tp. For the seven 
(9dsn > ^0 cases that were checked, the ratio Ar^i (square 
beam) : Ar» n jt (circular beam) was in the range 1.06-1.27, 
with a mean value of 1.15. 

B. Beam Offset Effects 

1. Results for Nonzero Integration Time. For a 
nonzero DSN-calibration beam offset [i.e., p / 0 in Eqs. (3) 
and (4)], the number of parameters (beamwidth, elevation 
angle, DSN-calibration offset distance and azimuth, and 
wind velocity and azimuth) becomes large. A simplified 
calculation was necessary in order to explore a reasonable 
volume of parameter space with the available computer 
resources. A pencil beam (zero diameter) was used for 
both the DSN antenna beam and the calibration beam. 
For this case, 6 ca \ — 9 dsn gives a delay difference with 
zero mean. 

For a given vector separation /f, the offset error depends 
on two azimuths: that of p and that of v w (in both cases, it 
is the azimuth relative to the DSN pointing direction which 
is important). Using pencil beams, the beam offset error 
was calculated for all combinations of these two azimuths 
(with grid spacings of 30 deg). Figures 3-5 show the results 
for separations (j|/7 1|) of 17, 50, and 100 m, and integration 
times of 30, 60, 120, and 240 sec. Figure 3 gi ves results 
for an elevation angle of 60 deg, Fig. 4 for an elevation 
angle of 20 deg, and Fig. 5 for an elevation angle of 6 deg. 
The lowest elevation angle currently used by the DSN for 
spacecraft tracking is 6 deg. The error bars in Figs. 3-5 
represent the standard deviation for the 144 azimuth pairs 
for each separation distance p and integration time. 


8 G. Lanyi, personal communication, Tracking Systems and Appli- 
cations Section, Jet Propulsion Laboratory, Pasadena, California, 
1992 . 
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For a 60-sec integration time and a maximum allowed 
error of 0.45 mm (i.e., a 20-percent contribution to a 1- 
mm 2 total variance), a separation of 50 m is acceptable at 
all elevation angles, and a separation of 100 m is acceptable 
at elevation angles larger than 15 deg. 

2. Tests of Numerical Results. In order to check 
the validity of the use of pencil beams for these calcu- 
lations, two cases were recalculated using more realistic 
beams (i.e., a 34-m-diam cylindrical DSN beam and a con- 
ical calibration beam). Both cases used 9 dsn — 10 deg, 
xp = 1 deg, an integration time of 100 sec, and a wind ve- 
locity of 8 m/sec in azimuth 45 deg. Case 1 had a 17-m 
DSN-calibration instrument separation in azimuth 90 deg, 
and case 2 had a 50-m separation in azimuth 0 deg. All 
azimuths are relative to the DSN pointing direction. The 
results agreed with the quadrature sum of the beam mis- 
match (i.e., the zero separation error) and the pencil-beam 
offset error to within 2 percent in both cases. 

The pencil beam results for a limited number of cases 
(16 combinations of elevation angle, separation distance 
and azimuth, wind azimuth, and integration time) were 
recalculated, using a different analytical derivation and 
computer code. In all 16 cases, the time- averaged beam 
offset error calculations agreed to <2 percent. 

IV. Results: Delay Rate Errors 

A. Beam Mismatch Effects 

The beam mismatch errors for delay rates [expressed as 
an Allan standard deviation: Eqs. (3) and (5) with /? = 0] 
were calculated for a cylindrical 34-m-diam DSN beam and 
a conical calibration beam as a function of elevation angle 
@dsn i calibration beam full width xp y and time interval At. 
Results are given in Table 1 (9 dsn = 60 deg) and Table 2 
(9 dsn = 20 deg). The values of xp = 1 deg and xp = 2 deg 
from Table 2 are plotted in Fig. 6. An elevation angle of 
20 deg is the lowest currently used for spacecraft Doppler 
tracking for navigation purposes, and is the lowest planned 
for use during the Cassini Gravitational Wave experiment. 
Column 6 in Tables 1 and 2 gives the Allan standard devi- 
ation of the uncalibrated wet troposphere for purposes of 
comparison. (It was calculated for a 34-m-diam cylindri- 
cal beam pointed in a constant direction at that elevation 
angle). The results in Tables 1 and 2 are averages over 
wind azimuth with a grid of 30 deg. The variation with 
wind direction was <5 percent for all cases except for the 
10-sec values at 20- deg elevation, where the la variation 
ranged from 5 percent at xp — 0.5 deg to 38 percent at xp = 
4 deg. In addition, the la variation for At — 32 sec, xj) — 
4 deg, and 9 dsn — 20 deg was 9 percent. 


For short time scales, other known error sources will 
dominate over the beam mismatch error for small beam- 
widths. Projected state-of-the-art frequency-standard ca- 
pability for the year 2000 9 is an Allan standard deviation 
of 1 x 10” 15 for At = 1-100 sec from a superconducting 
cavity maser oscillator and 1 x 10~ 14 /\/A t for At — 100— 
10,000 sec from a trapped-ion frequency standard. If water 
vapor radiometers are used for troposphere calibration (as 
seems likely), thermal receiver noise will be significant on 
short time scales. A WVR with a 200-MHz bandwidth 
and a 300-K total on-sky system temperature, operated in 
a total power mode, will give an Allan standard deviation 
of 7.4 x 10” 13 (A<) _l t“° tc 5 p , where tinteg is the WVR inte- 
gration time (in seconds) used for calibration. In Table 3, 
the quadrature sums of the Allan standard deviations for 
beam mismatch error, frequency standard noise, and WVR 
receiver noise (using the above parameters for the WVR 
and frequency standard) are given for an elevation angle 
of 30 deg and beamwidths of 0.5, 1, 2, and 4 deg. The 
assumed WVR integration time was tinteg = 10 sec for 
At < 100 sec, and tinteg = 100 sec for At > 100 sec. 
The improvement with decreasing beamwidth is substan- 
tial down to xp of approximately 1 deg for timescales 
up to At = 1000 sec. The improvement is weak for xp 
< 1 deg and At > 1000 sec. If future frequency stan- 
dards do not achieve the performance goals given above, 
delay-rate troposphere calibration will be limited by the 
frequency standards over a larger range of parameter space 
(e.g., for time scales shorter than 1000 sec and beamwidths 
wider than xp = 1 deg). 

B. Beam Offset Effects 

The delay-rate beam offset errors were calculated 
[Eqs. (3) and (5) with p / 0] for offsets of 17, 50, and 
100 m, using full beamwidths (xp) of 0.5, 1, 2, and 4 deg. 
The resulting Allan standard deviations were averaged 
over all wind directions (with a 30-deg grid) and over two 
DSN-calibration separation azimuths: 0 and 90 deg (as be- 
fore, all azimuths are relative to the DSN pointing direc- 
tion). Results (mean and standard deviation) are given in 
Tables 4 ( 9dsn = 60 deg) and 5 (9 dsn = 20 deg). The 17-, 
50- , and 100-m values from Table 5 are plotted in Fig. 7. 
Note that, for a 100-m offset, the beam offset error exceeds 
the uncalibrated troposphere for time scales less than ap- 
proximately the crossing time by the wind (100 m/8 m/sec 
= 12.5 sec). For longer time scales, the offset error de- 
creases, due to the transport of refractivity irregularities. 
Tables 4 and 5 list results for only one beamwidth: xp 
= 1 deg. Calculations for other beamwidths indicate a 


9 G. J. Dick, personal communication, Communications Systems Re- 
search Section, Jet Propulsion Laboratory, Pasadena, California, 
January 27, 1993. 
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weak dependence (<5 percent) upon beamwidth for sep- 
arations p > 50 m. For p = 17 m, the mean Allan stan- 
dard deviations varied by approximately ±10 percent with 
beamwidth at a 60-deg elevation angle. At p = 17 m and 
9 dsn = 20 deg, the mean Allan standard deviations for ip 
= 0.5, 2, and 4 deg are approximately 0.85, 1.5, and 2.5 
times that at ip = 1 deg. 

A DSN-calibration instrument offset has a severe im- 
pact on the accuracy of delay rate calibration. For a 50-m 
separation, which is probably near the minimum for a loca- 
tion off the antenna, the goals of the Cassini Gravitational 
Wave Search experiment cannot be met at time scales less 
than 3000-4000 sec, and the requirements cannot be met 
at time scales less than 500 sec. The use of multiple tro- 
posphere calibration instruments could potentially reduce 
the error. Due to time variation in wind azimuth, at least 
three instruments arrayed around the perimeter of a DSN 
antenna would be needed in order to achieve a significant 
reduction in the offset error. 

V. Discussion 

The consequences of a beam mismatch and offset be- 
tween a DSN antenna and a calibration instrument (e.g., a 
WVR) are much more serious for delay rate measurements 
than for time-averaged delay measurements. This situa- 
tion reflects the fact that differentiating a noisy signal (the 
wet tropospheric refractivity in this case) increases its am- 
plitude, whereas integrating such a signal decreases its am- 
plitude. For delay measurements with >l-min integration, 
HPBW <1.3 deg and an offset <50 m are needed for an 
overall calibration goal of 1 mm. A narrower beamwidth 
and an on-axis location would be desirable, in case other 
error sources should prove to be small enough to achieve 
a calibration accuracy of <1 mm. 

For delay rate measurements, the performance cost of 
an offset location for a calibration instrument is severe. A 


calibration instrument that could measure line-of-sight de- 
lay perfectly, but which had a 50-m offset from the axis of a 
DSN antenna, would calibrate only 96 percent of the total 
1000-sec troposphere Allan standard deviation, 80 percent 
for 100-sec intervals, and could do no useful calibration for 
10-sec time scales (these values are for a 20-deg elevation 
angle; the fractional calibration errors are approximately 
a factor of 2 smaller at 60 deg). An on- axis location for a 
calibration instrument is essential in order to perform ac- 
curate tropospheric delay rate measurements. Specifically, 
it is required in order to meet the Cassini Gravitational 
Wave Search goals at time scales less than 3000-4000 sec, 
or the requirements at times scales less than 500 sec. Multi- 
ple calibration instruments mounted around the perimeter 
of the DSN antenna could potentially reduce the offset er- 
ror, but this has not been quantified. Given the high cost 
of water- vapor-sensing instruments (e.g., WVR’s), this is 
not an attractive option. 

A calculation of the power spectral density error result- 
ing from an offset location reached conclusions similar to 
the results shown in Tables 4 and 5. 10 That power spec- 
tral density calculation was done in a very different manner 
from the calculations reported here (although it used the 
same wet troposphere fluctuation model [7]). 

The beam mismatch error (i.e., that for zero on-axis off- 
set) is important for time scales of 1000 sec and smaller. 
This error drops rapidly as the beamwidth is reduced, 
down to HPBW = 0.7 deg. The improvement for HPBW 
<0.7 deg is minor, due to contributions from other error 
sources (frequency standard noise and WVR thermal re- 
ceiver noise). For HPBW = 0.7 deg, the beam mismatch 
error will not interfere with the Cassini Gravitational Wave 
Search requirements at any time scale. It will prevent the 
goals from being met at time scales shorter than 200 sec. 


10 P. Richter, personal communication, Telecommunications Systems 
Section, Jet Propulsion Laboratory, Pasadena, California, 1993. 
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Table 1. Beam mismatch Allan standard deviation between a conical calibration beam and a 
coaxial 34-m-dIam DSN beam for a 60-deg elevation angle. 




Full beamwidth 


Uncalibrated 

Time interval 
(At), sec 


(V>), deg 


troposphere, 
Allan standard 
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1 2 

4 

deviation 
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X 
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X 
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X 
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X 
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1.1 

X 
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-16 

2.6 

X 
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1.6 

X 
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X 
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X 
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3.5 

X 
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-17 
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X 
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Table 2. Beam mismatch Allan standard deviation between a conical calibration beam and a 
coaxial 34-m-diam DSN beam for a 20-deg elevation angle. 




Full beamwidth 


Uncalibrated 

Time interval 
(At), sec 


(V0> deg 


troposphere, 
Allan standard 

0.5 

1 2 

4 

deviation 


10 

1.3 X 10~ 14 

32 

4.3 X 10" 15 
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8.3 X 10~ 14 

3.2 X 10~ 14 

7.5 X 10' 14 

1.1 X 10“ 14 

6.3 X 10“ 14 

3.5 X 10“ 15 

5.1 X IO -14 

1.1 X 10“ 15 

3.5 X 10“ 14 

3.5 X IO ” 16 

1.9 X 10~ 14 

1.1 X 10 “ lfi 

9.3 X 10“ 15 



TaW* 3. Total Allan standard deviation due to beam mismatch, frequency 
standard noise, and WVR receiver noise for a 30-deg elevation angle and an 
on-axis location. 
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Table 4. Allan standard deviation for an offset between a conical calibration beam and a 
34-m-diam DSN beam for a 60-deg elevation angle and a 1-deg full beam width. 


Time interval 

(At), sec _ 


DSN- calibration 
separation, m 


Uncalibrated 

troposphere, 

_ Allan standard 


17 

50 

100 

deviation 


10 

3.0 

± 0.4 

X 

10 " 
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X 
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X 
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32 
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-14 

4.6 
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X 
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-15 
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X 
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X 
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X 
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X 
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X 
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-16 
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X 
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10 - 
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X 
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X 
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X 

10 " 
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1.0 
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-16 
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10 - 

-16 

3.8 
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Table 5. Allan standard deviation for an offset between a conical calibration beam and a 
34-m-diam DSN beam for a 20-deg elevation angle and a 1-deg hill beam width. 


Time interval 
(At), sec 


DSN- calibration 
separation, m 


Uncalibrated 
trop osphere, 
Allan standard 
deviation 

17 

50 

100 

10 

3.4 ±1.8 x 10 -14 

7.1 ±4.1 x 10“ 14 

9.4 ±4.4 X 10 -14 

8.3 x 10“ 14 

32 

1.3 ±0.6 x 10~ 14 

3.6 ±1.5 x 10“ 14 

5.2 ±2.3 x 10- 14 

7.5 x lO" 14 

100 

4.8 ±2.1 x 10“ 15 

1.2 ±0.5 x 10“ 14 

2.0 ±0.8 x 10“ 14 

6.3 x 10“ 14 

320 

1.6 ±0.6 X 10- ,s 

4 .0 ±1.6 X 10~ 15 

7.2 ±2.6 X 10~ 1S 

5.1 X lO" 14 

1000 

5.2 ±2.0 X 10~ 16 

1.3 ±0.5 X 10“ 1S 

2.4 ±0.8 X 10~ 15 

3.5 X 10“ 14 

3200 

1.6 ±0.6 X 10~ 16 

4.0 ±1.6 X 10~ 16 

7.2 ±2.6 x 10- 16 

1.9 x 10~ 14 

10,000 

5.2 ±2.0 x 10~ 17 

1.3 ±0.5 X 10~ 16 

2.4 ±0.8 X 10- 1 ® 

9.3 X 10“ 15 
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Fig. 1. Instantaneous delay calibration error due to beam mis- 
match as a function of elevation angle ( 6 ) and full beamwidth (t/?) 
(circular beamwidth with a uniform gain) for a pencil DSN beam. 



1 10 100 
INTEGRATION TIME, sec 


Fig. 2. Time-averaged delay calibration error due to beam mis- 
match for four combinations of elevation angle (0) and full beam- 
width (V>) (circular beam with a uniform gain). A 34-m diam DSN 
beam was used in the calculations. 



Fig. 3. Time-averaged delay calibration error due to beam offset 
as a function of time for offsets (between the axis of the DSN 
antenna and the calibration Instrument) of 17, 50, and 100 m. The 
elevation angle Is 60 deg. A pencil DSN beam was used In the 
calculations. 



Fig. 4. Time-averaged delay calibration error due to beam offset 
as a function of time for offsets (between the axis of the DSN 
antenna and the calibration instrument) of 17, 50, and 100 m. The 
elevation angle Is 20 deg. A pencil DSN beam was used in the 
calculations. 
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Fig. 5. Time-averaged delay calibration error due to beam offset 
as a function of time for offsets (between the axis of the DSN 
antenna and the calibration instrument) of 17, 50, and 100 m. The 
elevation angle is 6 deg. A pencil DSN beam was used in the 
calculations. 



Fig. 6. Delay rate calibration error (expressed as an Allan stan- 
dard deviation) due to beam mismatch as a function of time In- 
terval for full beamwidths i/> (circular beam with a uniform gain) 
of V> = 2 deg and i/> = 1 deg. The elevation angle is 20 deg. 
The Allan standard deviation of the uncalibrated troposphere is 
plotted for comparison, as are the goals and requirements of the 
Cassini Gravitational Wave Search experiment. A 34-m-dlameter 
DSN beam was used In the calculations. 



Fig. 7. Delay rate calibration error (expressed as an Allan standard 
deviation) due to beam offset as a function of time interval for 
offsets (between the axis of the DSN antenna and the calibration 
instrument) of 17, 50, and 100 m. The elevation angle Is 20 deg. 
The Allan standard deviation of the uncalibrated troposphere is 
plotted for comparison, as are the goals and requirements of the 
Cassini Gravitational Wave Search experiment. A 34-m-dlameter 
DSN beam was used In the calculations. 
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